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ABSTRACT 

We present a strong and weak lensing reconstruction of the massive cluster Abell 2744, the first 
cluster for which deep Hubble Frontier Field (HFF) images and spectroscopy from the Grism Lens- 
Amplified Survey from Space (GLASS) are available. By performing a targeted search for emission 
lines in multiply imaged sources using the GLASS spectra, we obtain 5 high-confidence spectroscopic 
redshifts and 2 tentative ones. We confirm 1 strongly lensed system by detecting the same emission 
lines in all 3 multiple images. We also search for additional line emitters blindly and use the full 
GLASS spectroscopic catalog to test reliability of photometric redshifts for faint line emitters. We 
see a reasonable agreement between our photometric and spectroscopic redshift measurements, when 
including nebular emission in photometric redshift estimations. We introduce a stringent procedure to 
identify only secure multiple image sets based on colors, morphology, and spectroscopy. By combining 
7 multiple image systems with secure spectroscopic redshifts (at 5 distinct redshift planes) with 18 
multiple image systems with secure photometric redshifts, we reconstruct the gravitational potential 
of the cluster pixellated on an adaptive grid, using a total of 72 images. The resulting mass map is 
compared with a stellar mass map obtained from the deep Spitzer Frontier Fields data to study the 
relative distribution of stars and dark matter in the cluster. We find that the stellar to total mass ratio 
varies substantially across the cluster field, suggesting that stars do not trace exactly the total mass 
in this interacting system. The maps of convergence, shear, and magnification are made available in 
the standard HFF format. 

Subject headings: galaxies: evolution — galaxies: high-redshift — galaxies: clusters: individual (Abell 
2744) 


1. INTRODUCTION 

In the past two decades, gravitational lensing by clus¬ 
ters of galaxies has transitioned from an exotic curiosity 
to an invaluable tool for astrophysics and cosmology (e.g. 
Kneib & Natarajan 2011). Clusters can act as natural 
telescopes, magnifying background sources so that they 
appear brighter and more extended to the observer (e.g. 
Yee et al. 1996; Pettini et al. 2002; Bradac et al. 2012; 
Shirazi et al. 2014; Bayliss et al. 2014). The gravita¬ 
tional lensing effect can also be used to reconstruct the 
spatial distribution of mass in the clusters themselves, 
thus shedding light on the physics of dark matter and 
structure formation (e.g. Clowe et al. 2006; Bradac et al. 
2006; Sand et al. 2008; Newman et al. 2013; Sharon et al. 
2014; Merten et al. 2015). 

In the past two years, clusters of galaxies as tools 
for cosmology and astrophysics have become a major 
focus of a Hubble Space Telescope ( HST ) initiative. 
As part of the Hubble Frontier Field (HFF) program, 
six clusters of galaxies and six parallel fields are be¬ 
ing imaged to unprecedented depths in seven optical 
and near-infrared (NIR) bands, using the Wide Field 
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Camera 3 (WFC3) and the Advanced Camera for Sur¬ 
veys (ACS) (Coe et al. 2015). Similarly to previous 
public campaigns in legacy fields such as the Hub¬ 
ble Deep Fields (Williams et al. 1996; Ferguson et al. 
2000), this major HST campaign has triggered co¬ 
ordinated observations with all major facilities, rang¬ 
ing from Director’s Discretionary time on the Spitzer 
Space Telescope to deep observations with the Chan¬ 
dra X-ray Telescope and many ground based facilities 
(see the HFF website at http://www.stsci.edu/hst/ 
campaigns/frontier-fields/FF-Data for details). 

One of the major drivers of the HFF initiative is to 
search for magnified objects behind galaxy clusters. Ac¬ 
curate magnification maps are needed to determine the 
unlensed (intrinsic) properties of these galaxies. Fortu¬ 
nately, the deep images increase significantly the number 
of known multiply imaged systems per cluster, thus im¬ 
proving dramatically the fidelity of the mass models, and 
magnification maps. Typical CLASH 1 -like depth imag¬ 
ing of clusters (limiting magnitude ^27 ABmag for a 
5 -a point source Postman et al. 2012) yields a handful of 
multiple image systems per cluster (e.g. Johnson et al. 
2014), however at the depth of the HFF we expected 

1 http://www.stsci.edu/~postman/CLASH/Home.html 
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to find tens of multiply image systems. This was beau¬ 
tifully confirmed by the first papers that analyzed the 
HFF imaging datasets (e.g. Jauzac et al. 2014; Ishigaki 
et al. 2015; Laporte et al. 2015). 

The spectacular increase in the quality of imaging data 
has not yet been matched by advances in spectroscopic 
redshift (z spec ) determination for the multiple images. 
Thus, modelers often have to rely on photometric red- 
shifts to incorporate multiple images in their analysis, 
or sometimes they decide to leave the source redshift as 
a free parameter to be inferred by the lens model it¬ 
self. Whereas one of these two choices is inevitable in 
the absence of other data, it is also fraught with peril. 
Photometric redshifts can be very uncertain and prone to 
catastrophic errors for sources that are well beyond the 
spectroscopic limit, such as most of the faint arcs and 
arclets. Similarly, letting the mass model determine the 
redshift of the multiple images can potentially introduce 
confirmation bias in the modeling process. 

Obtaining spectroscopic redshifts for as many multiple 
image systems as possible is thus a fundamental step if 
we want to improve the mass models of the HFF clusters, 
and thus make the most of this groundbreaking initiative. 
Several efforts are currently underway to secure these 
redshifts, including our own based on the Grism Lens 
Amplified Survey from Space (GLASS) data, which is 
presented in this paper. GLASS is a large HST program 
that has just completed obtaining deep spectroscopy in 
the fields of ten clusters, including all six in the HFF 
program. A full description of the survey and its scientific 
drivers is given in paper I (Treu et al. 2015, submitted). 

In this paper we present a new mass model for the 
galaxy cluster Abell 2744, the first cluster for which both 
GLASS and HFF complete datasets are available. By ex¬ 
ploiting the exquisite imaging and spectroscopic data we 
carry out a rigorous selection of multiple images used 
to constrain the mass model. We use spectroscopic red¬ 
shifts - including 55 new line emitters detected at high- 
confidence - when available to supplement the photo¬ 
metric redshift (z p hot) calibration. We then compare the 
inferred two dimensional mass distribution to the distri¬ 
bution of stellar mass as determined from the analysis of 
deep Spitzer imaging data, indicating that the stellar to 
mass ratio varies substantially across the cluster, which 
is expected given the merging nature of this cluster. 

The paper is organized as follows. In Section 2 we give 
an overview of the data. In Section 3 we describe the 
reduction and analysis of the GLASS data. In Section 4 
we detail our algorithm for selection of multiple image 
systems. In Section 5 we present our mass model and 
study the relative distribution of stellar and total mass. 
In Section 7 we summarize our conclusions. Throughout 
this paper we adopt a standard concordance cosmology 
with Tt m = 0.3, = 0.7 and h = 0.7. All magnitudes 

are given in the AB system (Oke 1974). 

2. DATA 

Being one of the most extensively studied galaxy clus¬ 
ters, Abell 2744 is observed on a spectrum extending 
from X-ray to radio bands (e.g. Kempner & David 2004). 
In this paper, we focus on the optical and NIR imaging 
and spectroscopy data newly acquired with the Hubble 
and Spitzer Space Telescopes, as part of the GLASS pro¬ 
gram (2.1) and HFF initiative (2.2 and 2.3). 


2.1. The Grism Lens-Amplified Survey from Space 

The Grism Lens-Amplified Survey from Space 2 
(GLASS, GO-13459, P.I. Treu) is observing the cores 
of 10 massive galaxy clusters with the HST WFC3 NIR 
grisms G102 and G141 providing an uninterrupted wave¬ 
length coverage from 0.8/im to 1.7/im. The slitless spec¬ 
troscopy is distributed over 140 orbits of HST time in 
cycle 21. The last GLASS observations were taken in 
January 2015. Amongst the 10 GLASS clusters, 6 are 
targeted by the HFF (see Section 2.2) and 8 by the Clus¬ 
ter Lensing And Supernova survey with Hubble (CLASH; 
P.I. Postman, Postman et al. 2012). Prior to each grism 
exposure, imaging through either F105W or F140W is 
obtained to assist the extraction of the spectra and the 
modeling of contamination from nearby objects on the 
sky. The total exposure time per cluster is 10 orbits in 
G102 (with either F105W or F140W) and 4 in G141 with 
F140W. Each cluster is observed at two position angles 
(P.A.s) approximately 90 degrees apart to facilitate de¬ 
blending and extraction of the spectra. 

In concert with the NIR observations on the cluster 
cores the two parallel fields corresponding to the two 
P.A.s are observed with ACS’s F814W filter and G800L 
grism. Each parallel field has a total exposure time of 7 
orbits. These observations map the cluster infall regions. 

A key focus point of GLASS is the advancement and 
improvement of the lens models of the 10 surveyed clus¬ 
ters. This paper focuses on the modeling of the first 
cluster in the GLASS survey with complete GLASS spec¬ 
troscopy and HFF imaging as described in Section 3. 

2.2. Hubble Frontier Fields 

The Hubble Frontier Fields initiative 3 (HFF, P.I. Lotz) 
is a Director’s Discretionary Time legacy program with 
HST devoting 840 orbits of HST time to acquire opti¬ 
cal ACS and NIR WFC3 imaging of six of the strongest 
lensing galaxy clusters on the sky; Abell 370, Abell 2744, 
MACSJ 2129, MACSJ 0416, MACSJ 0717, and MACSJ 
1149. For a 5-cr point source, the limiting magnitudes 
are roughly 29 ABmag in both the ACS/optical (F435W, 
F606W, F814W) and WFC3/IR filters (F105W, F125W, 
F140W, F160W). All six HFF clusters are included in 
the GLASS sample described in Section 2.1. The pro¬ 
gram was initiated in HST cycle 21 and is planned for 
completion in 2016. The first cluster to have complete 
HFF and GLASS data available is the cluster studied in 
this paper, namely Abell 2744. 

An important aspect of the HFF efforts has been the 
community efforts to model the lensing clusters. Prior 
to the start of the HFF observations five independent 
groups (CATS team, Sharon et al., Merten, Zitrin et al., 
Williams et al., and Bradac et al.) provided lens mod¬ 
els of the HFF clusters, which have been made available 
online 4 . In the remainder of this work we will use our 
own models from Bradac et al., but we will make com¬ 
parison with several of the other publicly available HFF 
lens models of Abell 2744. 

2.3. Spitzer Frontier Fields 

2 http://glass.physics.ucsb.edu 

3 http://www.stsci.edu/hst/campaigns/frontier-fields 

4 http://archive.stsci.edu/prepds/frontier/lensmodels/ 
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The Spitzer Frontier Fields program 5 (P.I. Soifer) is a 
Director’s Discretionary Time program that images all 
six strong lensing galaxy clusters targeted by the HFF in 
both warm IRAC channels (3.6 and 4.5 /im). With the 
addition of archival imaging from IRA C Lensing Survey 
(P.I. Egami) we reduced the data using the methodol¬ 
ogy employed by the Spitzer UltRa Faint SUrvey Pro¬ 
gram (SURFS UP, P.I. Bradac, Bradac et al. 2014). 
IRAC imaging reaches ^ 50 hr depth on the primary 
(Abell 2744 cluster) and parallel (^6' to the west) fields. 
There are four flanking fields with ~ 25 hr depth (two 
to the north and two to the south) of the primary and 
parallel fields. The nominal depth for a 5-cr point source 
can reach 26.6 ABmag at 3.6 /am and 26.0 ABmag at 4.5 
/im, respectively. However this sensitivity might be com¬ 
promised in cluster center due to blending with cluster 
members and the diffuse intra-cluster light (ICL). 

For this work, we focus on the primary field of 
Abell 2744 that we have processed in a fashion very 
similar to that discussed in Bradac et al. (2014); Ryan 
et al. (2014); Huang et al. (2015). In brief, we applied 
additional warm-mission column pulldown and auto- 
muxstripe corrections as provided by the Spitzer Science 
Center to the corrected basic-calibrated data (cBCD) to 
improve image quality, particularly near bright stars. We 
process these files through a standard overlap correction 
to equalize the sky backgrounds of the individual cBCDs. 
We drizzle these sky-corrected cBCDs to 076 pix -1 with 
a DRIZ_FAC= 0.01 using the standard MOPEX soft¬ 
ware 6 . As a final note, there are ^ 1800 frames (with 
FRAMETIME^ 100 s) per output pixel in the deep regions. 

3. GLASS OBSERVATION AND DATA REDUCTION 

The two P.A.s of GLASS data analyzed in this study 
were taken on August 22 and 23 2014 (P.A. = 135) and 
October 24 and 25 2014 (P.A. = 233), respectively. Prior 
to reducing the complete GLASS data each exposure 
was checked for elevated background from He Earth-glow 
(Brammer et al. 2014) and removed, if necessary. The 
Abell 2744 data was favorably scheduled so only the Au¬ 
gust 23 reads were affected and thus removed by our 
reduction pipeline. 

The resulting total exposure times for the individ¬ 
ual grism observations are: G102_PA135 10929 sec¬ 
onds, G141.PA135 4212 seconds, G102.PA233 10929 sec¬ 
onds, and G141_PA233 4312 seconds. The correspond¬ 
ing exposure times for the direct GLASS imaging are: 
F105W_PA135 1068 seconds, F140W_PA135 1423 sec¬ 
onds, F105W.PA233 1068 seconds, and F140W.PA233 
1423 seconds. 

In Figure 1 we show a color composite image of 
Abell 2744, using the optical and NIR imaging from the 
HFF combined with the NIR imaging from GLASS. The 
red and green squares show the positions of the spectro¬ 
scopic GLASS Abell 2744 field-of-views. 

The GLASS observations are designed to follow the 
3D-HST observing strategy (Brammer et al. 2012). The 
data were processed with an updated version of the 3D- 
HST reduction pipeline 7 . Below we summarize the main 

5 http://ssc.spitzer.caltech.edu/warmmission/scheduling/ 
approvedprograms/ddt/frontier/ 

6 http://irsa.ipac.caltech.edu/data/SPITZER/docs/ 
dataanalysistools/tools/mopex/ 

7 http://code.google.com/p/threedhst/ 


steps in the reduction process of the GLASS data but 
refer to Brammer et al. (2012) and the GLASS survey 
paper (Treu et al. 2015) for further details. 

The data were taken in a 4-point dither pattern identi¬ 
cal to the one shown in Figure 3 of Brammer et al. (2012) 
to optimize rejection of bad pixels and cosmic rays and 
improve sampling of the WFC3 point spread function. 
At each dither position, a combination, a direct (F105W 
or F140W), and a grism (G102 or G141) exposure were 
taken. The direct imaging is commonly taken in the fil¬ 
ter “matching” the grism, i.e. pairs of F105W-G102 and 
F140W-G141. However, to accommodate searches for 
supernovae and the characterization of their curves in 
GLASS clusters, each individual visit is designed to have 
imaging in both filters. Hence several pairs of F140W- 
G102 observations exist in the GLASS data. This does 
not affect the reduction and extraction of the individual 
GLASS spectra. 

The individual exposures were turned into mosaics us¬ 
ing AstroDrizzle from the DrizzlePac (Gonzaga & et 
al. 2012), the replacement for MultiDrizzle (Koekemoer 
et al. 2003) used in earlier versions of the 3D-HST re¬ 
duction pipeline. Then all exposures and visits were 
aligned using tweakreg, with backgrounds subtracted 
from the direct images by fitting a second order polyno¬ 
mial to each of the source-subtracted exposures. We sub¬ 
tracted the backgrounds using the master sky presented 
by Kiimmel et al. (2011) for the G102 grism, while for 
the G141 data the master backgrounds were developed 
by Brammer et al. (2012) for 3D-HST. The individual 
sky-subtracted exposures were combined using a pixel 
scale of 0706 per pixel as described by Brammer et al. 
(2013) (^half a native WFC3 pixel) which corresponds to 
roughly 12A/pixel and 22A/pixel for the G102 and G141 
grism dispersions, respectively. Figure 2 shows these full 
field-of-view mosaics of the two NIR grisms (G102 on 
the left and G141 on the right) at the two GLASS P.A. 
for Abell 2744. The individual spectra were extracted 
from these mosaics by predicting the position and ex¬ 
tent of each two-dimensional spectrum according to the 
SExtractor (Bertin & Arnouts 1996) segmentation maps 
from the corresponding direct images. This is done for 
every single object and contaminations, i.e., dispersed 
light from neighboring objects in the direct image, so 
these contaminations can be estimated and accounted 
for. 

4. IDENTIFICATION OF MULTIPLE IMAGES 

In this section we describe how we identify and vet 
multiple image candidates using imaging (4.1) and spec¬ 
troscopic (4.2 and 4.3) data. 

4.1. Imaging data: identification and photometric 
redshifts 

Merten et al. (2011) published the first detailed strong 
lensing analysis of Abell 2744 identifying a total of 34 
multiple images (11 source galaxies) in imaging data 
pre-dating the HFF. With the addition of the much 
deeper HFF data, now a total of 176/56 candidate arc 
images/systems have been identified prior to this work 
(Atek et al. 2014, 2015; Zitrin et al. 2014; Richard et al. 
2014; Jauzac et al. 2014; Lam et al. 2014; Ishigaki et al. 
2015, see also Figure 3 and Table 1). We identify a new 
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Fig. l. — : The color composite image of Abell 2744 based on the HFF and GLASS imaging. The blue, green and 
red channels are composed by the filters on the right. The two distinct P.A.s of the spectroscopic GLASS pointings 
are shown by the red (P.A.=233 degrees) and green (P.A.=135 degrees) squares. Figure 2 shows the grism images 
corresponding to these two P.A.s. The locations of the emission line objects from Table 2 are marked by circles, with 
color coding reflecting the GLASS spectroscopic redshift quality (cf. Section 4.2 and column “Quality” in Table 2; 
2=blue, 3=orange, and 4=magenta). The cyan square shows the outline of Figure 3. 


multiply imaged system (i.e. system 60 in Table 1) com¬ 
prised of three images. 

Despite the vast number of identifications, only a hand¬ 
ful of multiple images have been spectroscopically con¬ 
firmed; arcs 4.3 and 6.1 were spectroscopically confirmed 
by Richard et al. (2014), arcs 3.1, 3.2 and 4.5 by John¬ 
son et al. (2014), and arc 18.3 by Clement et al. (in 
preparation). All corresponding redshifts can be found in 
Table 1. When spectroscopy is lacking, confirming that 
images belong to the same source is much more difficult. 
Photometric redshifts alone are not adequate for confir¬ 
mation. Ilbert et al. (2006) found that, within a given 
survey, the fraction of catastrophic errors in photometric 
redshift increases with faintness and redshift. Multiple 
images are typically faint and are necessarily at redshifts 
larger than the cluster redshift, which is z = 0.308 for 
Abell 2744. The mean F140W magnitude of all images in 
Table 1 is 27.11, and the mean source redshift is z = 2.63. 
Photometry of sources in cluster fields is complicated due 


to blending with cluster members and the ICL. While 
we do compute photometric redshifts in order to use the 
multiple images in our lens model, we do not rely on 
them alone to test the fidelity of the images. 

We instead rely on color and morphology information 
to determine whether images belong to the same mul¬ 
tiple image system; to first order, multiple images of 
the same source have identical colors. To compute col¬ 
ors, photometry is done using SExtractor in dual-image 
mode. We use F160W as the detection image because it 
detects the largest fraction of all multiple image candi¬ 
dates. We then measure isophot al magnitudes and errors 
in all seven photometric filters. Due to difficulties in de¬ 
tecting many of the multiple image candidates using the 
default SExtractor settings, we adopt a more aggressive 
set of settings for the objects with low S/N and/or highly 
blended. We refer to the default SExtractor settings as 
“cold” mode and the more aggressive one as “hot” mode 
photometry. These are similar in spirit but not identical 
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Fig. 2.—: The GLASS G102 (left) and G141 (right) grism pointings of Abell 2744 at two distinct P.A.s, with field- 
of-view shown by the red (P.A.=233 degrees) and green (P.A.=135 degrees) squares. The circles in all panels denote 
the positions of the emission line objects identified in this work, with color coding following the convention adopted 
in Figure 1. 


to those adopted by Guo et al. (2013). Even with the 
“hot” mode settings, we cannot detect all of the mul¬ 
tiple image candidates, though the detected fraction is 
vastly increased over the “cold” mode settings. Using 
the seven HFF photometric filters, we compute 4 colors 
for each image: F435W — F6061U, F814W — FT05IU, 
F125LF - F140VF, and FU0W - F160W 8 . The colors 
are computed within a fixed aperture (MAG_APER) that 
is 0.4" in diameter. We compute a reduced “color-y 2 ” 
for each image: 


xL = 



where i runs over the number of colors, N is the total 
number of colors we are able to measure, C{ is the i-th 
color, Ci is the inverse variance-weighted mean color of 


8 Note that the last two colors are not independent due to the 
repetition of F1AQW. We chose to repeat one filter to increase the 
number of color bins. 


all images in the system, and cq is the uncertainty in the 
i-th color. 

Multiple images of the same source also have pre¬ 
dictable morphologies. In rare cases, more than one im¬ 
ages of the same source possess a number of uniquely 
identifiable features. For instance, there are two such 
systems in Abell 2744, i.e., systems 1 and 2. The coun¬ 
terparts to systems 1 and 2 are systems 55 and 56, re¬ 
spectively. We include both counterpart systems in the 
lens model at the spectroscopic redshifts of systems 1 and 
2 that we measure in this work. All multiple images in 
Table 1 are visually inspected. We assign each image a 
grade that determines the likelihood that it is part of the 
system to which it is assigned. We perform this grading 
exercise in a lens model-independent fashion; we do not 
make any assumptions about the location of the critical 
curves relative to the graded images. There are, however, 
configurations of multiple images that are impossible to 
achieve through gravitational lensing by galaxy clusters, 
such as three individual images (not part of an elongated 
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Fig. 3.— : All multiple images discovered to date in Abell 2744. Green images are used in the lens model, while red 
images are unused for reasons discussed in Section 4. Cyan images belong to a new multiple image system discovered in 
this work. The yellow line is the critical curve from our best-fit lens model at z = 9. The RGB image is a combination 
of F105W, F125W and F160W. 


arc) on the same side of, and very distant from the clus¬ 
ter core with no counter-images. Other information such 
as surface brightness and symmetry can be incorporated 
independently of an assumed mass distribution, and we 
rely on this information much more heavily in assigning 
the morphology grade. The grading scheme based on 
morphological similarity is as follows: 

4 Image is definitely part of the system 

3 Image is very likely part of the system 

2 Image is potentially part of the system 

1 Image is very unlikely part of the system 

0 Image is definitely not part of the system 

Two inspectors (A.H. and M.B.) independently assign 
a grade to each image. The inspectors use several RGB 
images of the full HFF depth that span the full HST spec¬ 
tral coverage to assign the grade for each image. The two 
grades are then summed to get the reported morphology 
grade. Examples of multiple images that receive high 
and low morphology grades are shown in the Appendix. 

We use the color and morphology information together 
to determine whether to include a multiple image in our 
model. The joint criteria are: 

(xL <1-5 V M > 5) A M > 0, (2) 


where M is the summed morphology grade from each in¬ 
spector, which ranges from 0-8. In cases where the con¬ 
tamination by foreground objects, clusters members or 
ICL is severe, we rely only on the morphology criterion, 
M > 5. The particular threshold value of 1.5 was 
chosen because it is in the typical range for a good re¬ 
duced chi-square test, and most of the images in spectro¬ 
scopically confirmed systems are below this value. M > 5 
is chosen because in the least confident case that obeys 
this, M = 6, the modelers either both think the image is 
very likely part of the system (M = 3) or one thinks the 
image is potentially part of the system (M = 2), while 
the other is sure of it (M = 4). 

For a multiple image system to constrain the lens 
model, we must estimate its redshift. Having photo¬ 
metric redshift measurements for multiple images of the 
same source can provide a tighter constraint than a single 
measurement. Individual photometric redshifts are com¬ 
puted using EAZY (Brammer et al. 2008). We then use 
a hierarchical Bayesian model to obtain a single redshift 
probability density function for each system (Dahlen 
et al. 2013). The mode of this probability density func¬ 
tion will be referred to as ^Bayes- Only non-contaminated 
objects contribute to calculating ^Bayes- We graphically 
outline the procedure for measuring and including photo¬ 
metric redshifts as inputs to our lens model in Figure 4. 
2/57 systems (36 and 52) are entirely contaminated, so 
we do not compute ^Bayes for those systems. 14/57 sys¬ 
tems have ^Bayes < ^Cluster, and thereby are not included 
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in lens modeling (a fraction of those have poorly con¬ 
strained posteriors, monotonically declining from 0; they 
are highly uncertain and considered unreliable; we label 
them by assigning ^Bayes = 0.01). 5/57 systems have a 
multi-modal or extremely broad Bayesian redshift dis¬ 
tribution. We similarly do not include these systems in 
the lens modeling. For systems where only one image 
passed the color/morphology criteria and that image is 
not contaminated, we report ZBayes i n Table 1, but we do 
not include these systems in the modeling. We show the 
posterior redshift distributions for some of these cases in 
the Appendix. System 18 consists of a spectroscopically 
confirmed image (18.3), but the system as a whole does 
not pass all criteria required to be considered a multiple 
image system. Our screening rules out all the above- 
mentioned systems and delivers a secure set of 25/57 
multiple arc systems that are used in lens modeling. 

4.2. Targeted GLASS spectroscopy 

GLASS spectroscopy was carefully examined for a to¬ 
tal of 179 multiply lensed arc candidates mostly seen on 
both P.A.s with the goal of measuring spectroscopic red- 
shifts. As described in GLASS paper I, each spectrum 
was visually inspected by multiple investigators (X.W. 
and K.B.S.) using the GLASS Graphic User Interfaces 
(GUIs) GLASS Inspection GUI (GiG) and GLASS In¬ 
spection GUI for redshifts (GiGz) The results were then 
combined and a preliminary list of arcs with emission 
lines was drafted. In the end, another round of double¬ 
check by re-running GiGz was also executed to make 
sure no potential emission lines were missed. Following 
GLASS procedure, a quality flag was given to the redshift 
measurement: Q=4 is secure; Q=3 is probable; Q=2 is 
possible; Q=1 is likely an artifact. As described in paper 
I, these quality criteria take into account the signal to 
noise ratio of the detection, the probability that the line 
is a contaminant, and the identification of the feature 
with a specific emission line. For example, Q=4 is given 
for spectra where multiple emission have been robustly 
detected; Q=3 is given for spectra where either a single 
strong emission line is robustly detected and the redshift 
identification is supported by the photometric redshift, or 
when more than one feature is marginally detected; Q=2 
is given for a single line detection of marginal quality. As 
shown in Table 2, new spectroscopic redshifts were ob¬ 
tained for 7 images in total, corresponding to 5 systems. 
Among them, 5 high-confidence (with quality flags 3 or 4) 
spectroscopic redshifts were measured for arcs 1.3, 6.1, 

6.2, 6.3, 56.1. The spectra of these objects are shown 
in Figures 5-9. In particular for arc 6.1, our spectro¬ 
scopic redshift matches that reported by Richard et al. 
(2014), and we provide the first spectroscopic confirma¬ 
tion that 6.2 and 6.3 are images of the same system. We 
note that our measured redshift for arc 56.1 z spec = 1.20 
(Q=3; probable) differs signficantly from that given by 
Johnson et al. (2014) for arc 2.1 z spec = 2.2 (possible), 
even though the two systems are likely to be physically 
connected. Our measurement is based on three pieces 
of evidence. First, we detect a spectral feature in G141 
at both P.A.s (see Figure 9 for details) with sufficiently 
high signal to noise ratio to study its spectra shape. The 
feature is better described by a single line (identified by 
us as Ha at 2 = 1.2) rather than a doublet like [OIII] 
(Ay 2 = 2.4). Second, a line is marginally detected in one 


of the G102 spectra at exactly the wavelength expected 
for [OIII] z = 1.2. Third, the wide spectral coverage of 
our data and the data available in the literature rule out 
the possibility of the feature we see in G141 being other 
prominent lines such as Mgll and [Oil]. Taking into ac¬ 
count both the evidence and previous results, we assign 
a quality flag of Q=3 (probable). 

The uncertainty on our spectroscopic redshift measure¬ 
ments is limited by the resolution of approximately 50A 
and by uncertainties in the zero point of the wavelength 
calibration. By comparing multiple observations of the 
same object we estimate the uncertainty of our measure¬ 
ments to be in the order of Az ~ 0.01. This is sufficient 
for our purposes and we will not pursue more aggressive 
approaches to improve the overall redshift precision (e.g. 
Brammer et al. 2012). 

4.3. Blind search in GLASS data 

The targeted spectroscopy done in Section 4.2 could 
potentially miss some multiply imaged sources that are 
not identified photometrically. In order to increase the 
redshift completeness of emission line sources (both mul¬ 
tiply and singly imaged), we also conducted a blind 
search within the entire grism field-of-view. As a first 
step, one of us (T.T.) visually inspected all the 2D grism 
spectroscopic data, the contamination models, and resid¬ 
uals after contamination for each of the 2445 objects in 
the prime filed of Abell 2744 given by the GLASS cat¬ 
alog. This yielded a list of 133 candidate emission line 
systems that were later on inspected using the GLASS 
GUIs GIG and GiGz by two of us (X.W. and K.B.S.) 
to confirm emission line identifications and measure red¬ 
shifts. In order to search for previously undiscovered 
multiple images we inspected each set of objects with 
mutually consistent redshift. None of the sets of galaxies 
at the same spectroscopic redshifts are consistent with 
being multiply lensed images of the same source. Some 
of them are ruled out because of their position in the 
sky, while others are ruled out because their colors and 
morphology are inconsistent with the lensing hypothesis. 

Nonetheless, we compiled a list of singly imaged emis¬ 
sion line objects, consisting of 18, 16, and 34 quality 2, 3, 
and 4 spectroscopic redshift measurements respectively, 
which are color coded in Figure 1. Among them, the 
high-confidence (with quality flag 3 or 4; orange and ma¬ 
genta circles in Figure 1) emission line identifications are 
also included in Table 2. As mentioned in Section 4.1, via 
running EAZY on the full-depth seven-filter HFF imag¬ 
ing data, we were able to measure photometric redshifts 
for those objects as well. As a result, a comparison be¬ 
tween spectroscopic and photometric redshift measure¬ 
ments is possible, as displayed by Figure 10. We find 
that 25/55 photometric redshifts agree within their 1-cr 
uncertainties with corresponding spectroscopic redshifts, 
when nebular emission lines are included in the fitting 
template. This suggests the presence of additional sys¬ 
tematic errors that are likely related to the photometric 
redshift fitting method. In order to account for the un¬ 
known systematics, we increase the photometric redshift 
uncertainties for the sources used in the construction of 
the lens model. 

We double-checked our spectroscopic redshift measure¬ 
ments by re-running GiGz on all the objects and also 
cross-checked the photometric redshift measurements 
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Fig. 4.— : Flow chart describing our procedure for assigning photometric redshifts to multiple image systems. Two 
segmentation maps, a “hot” and a “cold” version, were used for source detection. The detection and deblending 
thresholds are set to the SExtractor defaults in the cold segmentation map. Objects detected in the cold segmentation 
map, typically the brighter and more isolated ones, have more accurate photometry. The hot segmentation map was 
created using extremely aggressive detection and deblending thresholds. It detects the majority of the remaining 
multiple images that are not detected in the cold version. C/M cuts refer to the color/morphology cuts used to purify 
the sample of multiple images. ^Bayes refers to the redshift obtained from combining multiple redshifts in a hierarchical 
Bayesian model (Dahlen et al. 2013). ^Bayes is considered reliable if it is larger than the redshift of Abell 2744, 
z = 0.308, and it is not multi-valued. See the Appendix for examples of multiple image systems that pass and fail 
some of the tests in this flow chart. 


through re-fitting the photometric redshifts using a dif¬ 
ferent method by a subset of the authors (R.A., M.C., 
E.M) without knowing previous results. The general con¬ 
clusions about the agreement between spectroscopic and 
photometric analyses remains unchanged. 

5. GRAVITATIONAL LENS MODEL 

Our lens modeling method, SWUnited (Bradac et al. 
2009, Bradac et al. 2005), constrains the gravitational 
potential within a galaxy cluster field via x 2 minimiza¬ 
tion. It takes as input a simple initial model for the 
potential. A x 2 is then calculated from strong and weak 
gravitational lensing data on an adaptive, pixelated grid 
over the potential established by the initial model. The 
number of grid points is increased and the x 2 is recalcu¬ 
lated. Once the minimum is found, and convergence is 


achieved, derivative lensing quantity maps, such as con¬ 
vergence (ft), shear ( 7 ) and magnification (/i), are pro¬ 
duced from the best-fit potential map. Errors in these 
quantities are obtained via the method described below. 
Maps of the convergence and magnification are shown in 
Figure 11. 

A previous model of Abell 2744 using pre-HFF data 
was created using the same lens modeling code. The 
model was created as part of a call by STScI to model 
the HFF clusters, and it appears on the publicly accessi¬ 
ble HFF lens modeling website as the Bradac vl models 9 . 
The previous model was constrained using 44 total mul¬ 
tiple images belonging to 11 distinct systems. The weak 

9 http://www.stsci.edu/hst/campaigns/frontier-fields/ 
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Fig. 5.—: Emission line detection results on object ID 7^00963 (arc 6.1) at the two P.A.s displayed in the two 
sub-figures accordingly. In each sub-figure, the two panels on the first row show the observed 1-dimensional spectra, 
where the contamination subtracted flux is denoted by the blue solid line and the noise level by the cyan shaded 
region. The two panels on the second row give the corresponding contamination model in red dashed line. For the 
four panels directly underneath, the top two display the interlaced 2-dimensional spectra whereas the bottom two have 
contamination subtracted. In the 1- and 2-dimensional spectra, the identified emission lines are denoted by vertical 
dashed lines in magenta and arrows in red, respectively. The wheat colored regions cover contamination model defects. 
The two panels on the far left refer to the 2-dimensional postage stamp created from the HFF co-adds through drizzling 
(on top) and the 1-dimensional collapsed image (on bottom). Note that these two panels share the same x-axis along 
the grism dispersion direction. Some ancillary information can also be seen in the upper left corner in each sub-figure. 
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Fig. 6.— : Same as Figure 5, except that object ID 7^00523 (arc 6.2) is shown. 
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Fig. 7.— : Same as Figure 5, except that object ID ^00433 (arc 6.3) is shown. Note here WFC3/G102 at the second 
P.A. is cut off on the right due to grism defect. 



















































































































normalized flux normalized flux 


12 


Wang et al. (2015) 


WFC3/G102 


WFC3/G141 



0.10 - 

Abell 2744 

0.05- 

PA = 135 

o.oo- 

ID #00336 

| -0.05 ■ 
1/5 

O) -0.10 | 

^SDec = 1-498 

0 


Quality = 3 m 



Qy V 

arcseconds 



0.95 1.0 

A [/xm] 


1.3 1.4 

A [fim] 


WFC3/G102 


WFC3/G141 


Abell 2744 


0.0 

PA = 233 


-0.2 


C^l 

-0.4 


E 

ID #00336 

.0 

1/5 

-0.6 



-0.8 

* spec =1-498 

0 

-1.0 

Quality = 3 

0 

-1.2 



8 

6 




0.95 1.0 

A [/xm] 


1.3 1.4 

A [/xm] 


Fig. 8.— : Same as Figure 5, except that object ID ^00336 (arc 1.3) is shown. Moreover, the 2-dimensional spectra 
are smoothed, whilst the 1-dimensional spectral resolution remains unchanged. 
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Fig. 9.— : Same as Figure 8, except that object ID 7^00888 (arc 56.1) is shown. 
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Fig. 10. — : Comparison between the spectroscopic and pho¬ 
tometric redshifts for the 55 objects with high-confidence emission 
lines (quality flags 3 or 4). We also show the 1-cr error bars (enclos¬ 
ing 68% of the total probability) around the photometric redshifts. 
There is reasonably good agreement between photometric and spec¬ 
troscopic redshifts, with 25 out of 55 spectroscopic redshifts within 
the photometric redshift error bars. This is acceptable consider¬ 
ing that the photometric redshift uncertainties only include the 
random component, even though they suggest that an additional 
systematic uncertainty component is present. In order to account 
for this systematic uncertainty, for the systems used to build the 
lens model we added 20% in quadrature to our photometric errors. 

lensing constraints were obtained by one of the mod¬ 
elers, Julian Merten, and distributed to all participat¬ 
ing modeling teams. The same weak lensing constraints 
were used in the model that appears in this work. This 
model is also made available to the public on the HFF 
lens modeling website as the Bradac v2 model. In the 
vl model, magnification uncertainties were estimated by 
bootstrap-resampling the weak lensing galaxies. In this 
work, however, we took a different approach to estimate 
uncertainties, one that we expect more accurately rep¬ 
resents the true uncertainties. Because the number of 
multiple image systems used in this model is much larger 
than in the vl model, 72 total multiple images belong¬ 
ing to 25 distinct systems, we bootstrap-resampled the 
multiple image systems that were not spectroscopically 
confirmed. These are the systems for which we use ^Bayes 
in the lens model. We assess the impact of photomet¬ 
ric redshift uncertainty on the derived lensing quanti¬ 
ties by resampling the redshift of each system lacking 
spectroscopic confirmation from their full ^Bayes poste¬ 
riors 10 . We compare the variance in magnification due 
to redshift uncertainty with the variance in magnification 
due to bootstrap-resampling the multiple image systems, 
finding that the latter is dominant. We nonetheless prop¬ 
agate both sources of error when reporting the errors on 
all derived lensing quantities in this work. 

As a test of the improvement of the lens model with 
the addition of the new multiple image constraints from 
the HFF data, we calculate the magnification of SN 
HFF14Tom, a Type la Supernova (SN la) at z = 1.35 dis¬ 
covered in the primary cluster field of Abell 2744 (Rod¬ 
ney et al. 2015). We compare the magnification predicted 

10 We exclude values of the redshift z < Cluster + 0-1 when 
resampling from the ^Bayes posteriors because they are unphysical. 


by our lens model with the magnification calculated di¬ 
rectly from a comparison with other SNe la at similar 
redshifts, /i = 2.03 =b 0.29 (Rodney et al. 2015). The 
magnification predicted by the vl model, using pre-HFF 
data was /u be st = 3.15, with /ti me dian = 2.45jjfie (68% 
confidence). The new model presented in this work, v2, 
predicts a consistent magnification of fi best = 2.23, with 
/^median = 2.24+o'°g (68% confidence). The improved 
lensing constraints significantly improve the accuracy as 
well as the precision according to this test. We note that 
while we were not blind to the magnification of the super¬ 
nova predicted by Rodney et al. (2015) when producing 
the v2 lens model, we did not use the magnification as 
an input to our model. 

5.1. Comparison with previous work 

Three teams (Jauzac et al. 2014, Lam et al. 2014 and 
Ishigaki et al. 2015) have published models of Abell 2744 
using new multiple image constraints identified in the 
HFF imaging data. Of these teams, currently only the 
lens models produced by Ishigaki et al. 2015 (GLAFIC) 
are publicly available through the Mikulski Archive for 
Space Telescopes (MAST 11 ). We compare our models to 
theirs as well as the Sharon v2 models, which include 
updated spectroscopic measurements of multiple images 
identified before the HFF data were obtained (Johnson 
et al. 2014). Finally, we also compare our model with 
several updates of the CATS(vl) models (Jauzac et al. 
2015, private communication). The CATS(v2) models 
are presented by Jauzac et al. (2014) and use a much 
larger number of multiple images than we include in our 
lens model. CATS(v2.1) employs the same set of multi¬ 
ple images as CATS(v2), but makes use of the spectro¬ 
scopic redshifts obtained in this work. CATS(v2.2) uses 
the same set of multiple image constraints used in our 
model. We compare the surface mass density profiles 
(Figure 12) and cumulative magnified source plane ar¬ 
eas (Figure 13) predicted by all models described above. 
The surface mass density profiles agree quite well at radii 
where multiple image constraints are plentiful. However, 
the models begin to differ rapidly near the boundaries 
of the constrained area. Our model disagrees with the 
CATS models most severely. It is interesting to note that 
the three CATS models agree internally extremely well, 
despite CATS(v2.2) using the same set of multiple images 
used in this work, a considerably different set than the 
one used in CATS(v2) and CATS(v2.1). The significant 
difference between our model and the CATS models may 
be due to differences in the modeling techniques or the 
fact that our method uses additional constraints (weak 
lensing). Weak lensing constraints have a stronger im¬ 
pact on the model at radii beyond where multiple images 
are observed. In contrast, there is excellent agreement 
among the models in the inferred magnified source plane 
area. Thus, even though there may be small but sig¬ 
nificant differences in the specific details of each recon¬ 
structions, by and large the total integrated properties 
are very similar. 

We also note that our model supersedes the model ob¬ 
tained by members of our team as part of the initial HFF 
modeling effort based on pre-HFF data. The uncertain¬ 
ties in this current version of the model are smaller, con- 

11 http://archive.stsci.edu/prepds/frontier/lensmodels/ 
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Fig. 11.—: Convergence (left), k, and flux magnification (right), /r, maps of Abell 2744 produced by our lens model for a source at z = 9. 
Maps cover 3.5 X 3.5arcmin 2 . In both maps, the gray stick pattern indicates the phase angle of the shear. The figures are made using 
software written by Dan Coe. 
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Fig. 12.— : Surface mass density profile for the lens 
model obtained in this work compared to several re¬ 
cently published lens models of Abell 2744. The Sharon 
v2 model is presented by Johnson et al. (2014), the 
GLAFIC model by Ishigaki et al. (2015), the CATS(v2) 
model in Jauzac et al. (2014), and the CATS(v2.1) and 
CATS(v2.2) models by Jauzac et al. (2015, private com¬ 
munication). The shaded gray region indicates the radii 
over which multiple image constraints are available. The 
models agree best within this region, and they begin to 
significantly disagree at radii > 200kpc. The radius is 
measured from the center of the BCG. Error bars shown 
for our model represent 68% confidence. Gaussian 1-cr 
error bars are included on all three CATS models, but 
are almost entirely too small to discern. 

sistent with the fact that we have increased the number 
of strong lensing constraints. 

We also compare our method of estimating redshifts of 
multiple image systems with the one used by the CATS 
team (Richard et al. 2014, Jauzac et al. 2014). In Fig¬ 
ure 14, ^Bayes is the redshift obtained from hierarchical 
Bayesian modeling of the photometric redshifts obtained 
in this work. z mo dei is the redshift obtained by Jauzac 
et al. (2014) by minimizing their analytical model un¬ 


Fig. 13.— : Cumulative source plane area versus magni¬ 
fication at z = 9. The models used in this comparison 
are the same as those described in Figure 12. 

certainty while leaving the redshift as a free parameter. 
It is important to check this procedure independently 
since leaving 2: mo dei as a free parameter or predicting ad¬ 
ditional multiple images that belong to the same system 
could in principle lead to confirmation bias. Overall, we 
find good agreement between ^Bayes and z m odeb within 
the admittedly large uncertainties on ^Bayes- There are 
only two new systems with spectroscopic redshifts avail¬ 
able to compare with z mo d e i, and they are both inconsis¬ 
tent at > 5 a. This may be due to small number statis¬ 
tics or perhaps an indication that the uncertainties on 
^modei are underestimated. More spectroscopic redshifts 
are needed to perform this test in a more stringent man¬ 
ner. 

6. THE SPATIAL DISTRIBUTION OF STELLAR AND DARK 

MATTER 

6.1. Stellar mass map 

The Spitzer IRAC 3.6 /am image samples close to rest- 
frame iGband for the cluster, so we use the 3.6 /am fluxes 
from cluster members to approximate the cluster stellar 
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^Bayes 


Fig. 14.— : Comparison of the redshifts determined in 
this work (^Bayes) versus the model-predicted redshifts 
given by Jauzac et al. (2014) for all multiple image sys¬ 
tems used in the lens model. Note that the previously 
confirmed spectroscopic systems are left out of this com¬ 
parison because z mo d e i were not calculated. Two systems 
(green) are spectroscopically confirmed in this work for 
the first time and are included in the comparison. For 
these two objects, we use the new spectroscopic redshift 
on the horizontal axis in place of ^Bayes- The 2: mo dei val¬ 
ues are in significant disagreement with the spectroscopic 
values for these two systems. ^Bayes represents the peak 
of a statistical combination of all available photomet¬ 
ric redshift probability density functions (Dahlen et al. 
2013). The vertical error bars reflect 1 -a Gaussian error 
on z m odei, and the horizontal error bars show the 68% 
credible interval for ^Bayes- 12/16 systems are consistent 
at 68% confidence level. 

mass distribution. We first selected the red sequence 
cluster members brighter than the 25th mag in F814W 
from the color-magnitude and color-color diagrams fol¬ 
lowing the procedure described in Richard et al. (2014). 
We also cross-matched the selected cluster members with 
the spectroscopic redshift catalog given by Owers et al. 
(2011) to ensure that we included all the cluster members 
confirmed with spectroscopy. We selected a total of 190 
bright cluster members for their stellar mass distribution. 

To create an image with 3.6 jam flux from cluster mem¬ 
bers only, we first created a mask with value 1 for pixels 
that belong to cluster members in the F160W image and 
0 otherwise. We then convolved the mask with the 3.6 fi m 
PSF to match the IRAC angular resolution, set the pixels 
below 10% of the peak value to zero, and resample the 
mask onto the IRAC pixel grid. We obtained the 3.6 /im 
map of cluster members by setting all IRAC pixels not 
belonging to cluster members to zero and smoothed the 
final surface brightness map with a two-pixel wide Gaus¬ 
sian kernel. 

The IRAC surface brightness map was transformed 


into a surface mass density map by transforming the 
3.6/im flux into K-band luminosity and then by multi¬ 
plying by stellar mass to light ratio derived by Bell & 
de Jong (2001) using the so-called “diet”-Salpeter stellar 
initial mass function (IMF). The resulting stellar mass 
map in show in the left panel of Figure 15. 

The main source of uncertainty on the stellar mass den¬ 
sity is the unknown initial mass function. For example, if 
one were to adopt a Salpeter (1955) IMF - as suggested 
by studies of massive early-type galaxies, the stellar mass 
density would increase by a factor of 1.55. 

6.2. Stellar to total mass ratio 

We obtain the stellar to total mass ratio map by divid¬ 
ing the stellar surface mass density map obtained from 
photometry by the total surface mass density map ob¬ 
tained from gravitational lensing. This is shown in the 
right panel of Figure 15. We note that resolution ef¬ 
fects are non-trivial to take into account, since the reso¬ 
lution of the lensing map depends on the density of local 
sources and the amount of regularization. Thus, the map 
should be interpreted keeping in mind this caveat. Inter¬ 
estingly the stellar to total mass ratio varies significantly 
across the cluster. Many but not all the massive ellip¬ 
ticals seem to reach values of 0.05 or more, which are 
typical of the central regions of isolated massive galaxies 
(e.g. Gavazzi et al. 2007). However, the ratio appears 
to be significantly lower in the center of the cluster and 
in the south-east quadrant. In future work, we plan to 
compare the observed map with those obtained from nu¬ 
merical simulations by taking into account the effects of 
finite resolution in the observed mass and light maps, in 
order to test whether the spread in stellar to total mass 
ratio is reproduced. Furthermore, we plan to carry out a 
systematic comparison with mass reconstructions where 
mass is assumed to follow light up to a scale factor (Zitrin 
et al. 2009). At face value our result is inconsistent with 
this assumption for a merging cluster like Abell 2744. 
However, uncertainties on both models and resolution 
effects must be taken into account in order to evaluate 
the significance of this apparent violation. Thus, this re¬ 
sult should be considered as preliminary until confirmed 
by a more detailed analysis. 


7. CONCLUSIONS 

In this paper we have used spectroscopic data from 
the GLASS survey in combination with ultra-deep imag¬ 
ing data from the HFF program to construct a strong 
gravitational lens model for the cluster Abell 2744. In 
an effort to obtain a precise and accurate mass model we 
carried out a systematic search for spectroscopic redshifts 
of multiple images and we applied a rigorous algorithm 
to select only secure multiple image systems amongst the 
dozens that have been proposed in the literature. The 
lensing mass map is then combined with a stellar mass 
map derived from IRAC photometry to study the rel¬ 
ative spatial distribution of luminous and dark matter. 
Our main results can be summarized as follows: 

1. We have measured spectroscopic redshifts for 5 
multiple image systems (quality flag 4 and 3, i.e. 
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Fig. 15.— : Stellar mass surface density (left) and stellar to total 
mass surface density is in the unit of M 0 kpc -2 . 

secure and probable). We have also confirmed spec¬ 
troscopically that images 6.1, 6.2, 6.3 belong to the 
same source. The spectroscopically confirmed im¬ 
ages are used to constrain the gravitational lens 
model. We also obtain 2 tentative redshifts, which 
are not used to to constrain the mass model, but 
could potentially be confirmed by future work. 

2. From the GLASS data we derive an extensive red- 
shift catalog of faint emission line systems which we 
use to test photometric and lensing determinations 
of redshift. Generally speaking, the measurements 
agree within the 1-a uncertainties, when nebular 
emission lines are included in the fitting template. 

In addition, we compare photometric redshifts with 
redshifts determined by Jauzac et al. (2014), based 
on their gravitational lens model and find an agree¬ 
ment within the large uncertainties of the former. 

For the two systems with new spectroscopic red¬ 
shifts we find a significant difference with respect to 
model redshifts. This may be due to small number 
statistics or to the model redshifts uncertainties be¬ 
ing underestimated. More spectroscopic redshifts 
are needed to make a more stringent test. 

3. Our rigorous selection algorithm identifies a to¬ 
tal of 25/72 multiple arc systems/images as secure 
out of a sample of 57/179 candidate multiple arc 
systems/images, compiled from the literature and 
from our own work. Most systems are rejected 
either on the basis of inconsistent morphology or 
inconsistent spectral energy distribution between 
the candidate multiple images, or because of in¬ 
sufficient evidence that they belong to the same 
source. 

4. The derived mass model is found to be very precise, 
as measured by bootstrap- and redshift-resampling 
the set of multiple images used as input. Further¬ 
more, we tested how well our model reproduces 
the magnification of the background SN la Tomas 
(Rodney et al. 2015). The SN la was not used as 
a constraint to the model and yet its magnifica¬ 
tion is consistent (2.03 db 0.29 for the supernova vs 
2 .231 q!o7 f rom our mass model). 


mass ratio (right) distributions of Abell 2744. Stellar 

5. Abell 2744 is confirmed to be an excellent gravita¬ 
tional telescope, with a source plane area of ~ 0.7 
arcminute square being magnified by a factor of 2. 

6. We construct a stellar surface mass density map 
and the stellar to total mass ratio by selecting the 
light associated with red sequence cluster galax¬ 
ies and using the total mass density map obtained 
from strong lensing. Albeit with significant uncer¬ 
tainties, we find that the stellar to mass ratio varies 
significantly across the cluster, tentatively suggest¬ 
ing that stellar mass does not trace total mass in 
this interacting system. 
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TABLE l:Multiply lensed arc systems identified in the Abell 2744 field 


IDarc 1 

R.A. 

Dec. 

F140W 

-^spec 

2 

^Bayes 



Quality 


Pass Flag 

Magnification 4 


(deg.) 

(deg.) 

(mag) 
















3 

^ Bayes 

Color-x 2 

Morphology 

Contaminated? 

System 

Image 


1.1 J 

3.597540 

-30.403920 

22.16 


1.70 [1.08-2.01] 


0.487 

8 

0 

1 

1 

3.88 [4.04-5.26] 

1.2 J 
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-30.406820 

22.24 




3.079 

8 

0 


1 
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1.3 G 
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1.50 
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8 

0 


1 

3.73 [3.52-4.66] 

2.1 J 

3.583250 

-30.403350 

24.87 

1.20 

0.01 


0.892 

8 

0 

1 

1 

7.34 [7.49-9.18] 

2.2 5 J 

3.597290 

-30.396720 






8 

0 


1 

2.53 [2.45-2.59] 

2.3 J 

3.585420 

-30.399900 

27.46 




1.702 

8 

1 


1 

14.84 [9.88-17.08] 

2.4 J 

3.586420 
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0.696 

8 

1 


1 
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3.1 T 
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-30.393870 

24.28 
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8 

0 

1 

1 
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8 

0 


1 
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0 


1 
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8 
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1 
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1 
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TABLE l:Multiply lensed 

arc 

systems identified 

in 

the Abell 2744 field 




25.1 J 

3.594460 

-30.402740 

27.23 

0.01 

0 

3.212 

6 

0 

0 

0 


25.2 J 

3.592170 

-30.403330 

27.38 



1.070 

6 

1 


0 


25.3 J 

3.584210 

-30.408290 

29.13 
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6 

0 


0 


26.1 J 
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28.28 

0.20 

0 
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5 

0 

0 

0 


26.2 J 
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0.809 

5 

1 


0 


26.3 J 
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5 

0 


0 


27.1 J 
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0.01 

0 
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5 

0 

0 

0 


27.2 J 
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5 

0 


0 


27.3 J 
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5 

1 


0 


28.1 J 
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27.20 

6.39 [3.68-7.44] 

1 

0.957 

7 

0 

1 

1 

4.88 [4.29-5.11] 

28.2 J 
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7 

0 


1 

3.14 [3.13-3.42] 

28.3 J 

3.585290 
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0.610 

6 

1 


1 

2.86 [2.89-4.13] 

28.4 J 

3.587460 
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0.181 

6 

1 


1 

3.24 [2.99-4.32] 

29.1 J 
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29.30 

4.82 [1.82-6.98] 

1 

0.394 

5 

0 

1 

1 

47.33 [19.69-103.88] 

29.2 J 

3.580540 

-30.400430 

28.97 



0.193 

5 

0 


1 

7.87 [6.10-8.44] 

29.3 J 

3.583630 

-30.396600 

30.58 



2.121 

4 

1 


0 

20.38 [21.25-50.53] 

30.1 J 

3.591000 

-30.397440 

25.77 

0.87 [0.75-6.47] 

1 

1.513 

8 

0 

1 

1 

8.21 [4.77-10.53] 

30.2 J 

3.586710 

-30.398190 

26.41 



3.772 

6 

1 


1 

4.52 [4.22-4.77] 

30.3 J 

3.581920 

-30.401700 

25.81 



0.914 

8 

0 


1 

4.07 [4.16-8.70] 

31.1 J 

3.585920 

-30.403170 

28.91 

11.90 

0 

0.507 

8 

1 

0 

0 


31.2 J 

3.583710 

-30.404120 

28.63 



0.398 

8 

1 


0 


31.3 J 

3.599830 

-30.395520 





3 

0 


0 


32.1 J 

3.583580 

-30.404720 

27.70 

0.71 [0.99-6.51] 

1 

0.769 

5 

0 

1 

1 

3.22 [2.71-3.30] 

32.2 J 

3.586670 

-30.403350 

27.43 



0.893 

5 

1 


0 

29.09 [16.59-26.21] 

32.3 J 

3.599790 

-30.395980 

26.78 



0.640 

4 

0 


1 

1.66 [1.65-1.71] 

33.1 J 

3.584710 

-30.403150 

28.72 

4.43 [3.12-5.54] 

1 

2.351 

7 

1 

1 

1 

667.25 [61.64-609.99] 

33.2 J 

3.584420 

-30.403390 

27.91 



0.356 

5 

0 


1 

16.15 [20.27-214.71] 

33.3 J 

3.600670 

-30.395420 

28.63 



0.664 

7 

0 


1 

2.39 [2.40-2.53] 

34.1 J 

3.593420 

-30.410840 

27.69 

3.82 [1.03-4.33] 

1 

1.898 

8 

0 

1 

1 

74.50 [9.39-37.41] 

34.2 J 

3.593830 

-30.410720 

27.84 



1.135 

5 

0 


1 

66.45 [11.48-40.75] 

34.3 J 

3.600580 

-30.404530 

26.91 



0.215 

6 

0 


1 

4.26 [3.98-5.10] 

35.1 J 

3.581080 

-30.400210 

27.00 

0.01 

0 

0.307 

7 

0 

0 

0 


35.2 J 

3.581540 

-30.399390 

26.33 



0.773 

7 

0 


0 


35.3 J 

3.597830 

-30.395540 

27.05 



2.534 

6 

0 


0 


36.1 J 

3.589460 

-30.394410 

29.50 

0.01 


4.101 

8 

1 

0 

0 


36.2 J 

3.588670 

-30.394300 

28.60 



2.210 

8 

1 


0 


36.3 J 

3.577500 

-30.401500 

28.28 



1.193 

6 

1 


0 


37.1 J 

3.589040 

-30.394910 

27.15 

0.01 


5.976 

7 

1 

0 

0 


37.2 J 

3.588710 

-30.394850 





7 

1 


0 


37.3 J 

3.578330 

-30.401400 

25.38 



2.109 

5 

0 


0 


38.1 J 

3.589420 

-30.394110 

29.12 

0.78 [0.99-5.63] 

0 

0.243 

6 

0 

0 

0 


38.2 J 

3.588960 

-30.394040 

28.39 



0.336 

6 

0 


0 


38.3 J 

3.576380 

-30.402130 

28.17 



1.254 

4 

0 


0 


39.1 J 

3.588790 

-30.392530 

27.61 

3.58 [0.80-4.10] 

1 

1.154 

7 

0 

1 

1 

12.03 [9.80-35.06] 

39.2 J 

3.588540 

-30.392510 

27.13 



1.355 

7 

0 


1 

53873.97 [21.17-147.82] 

39.3 J 

3.577460 

-30.399560 

25.49 



0.231 

2 

0 


1 

4.38 [3.79-5.90] 

40.1 J 

3.589080 

-30.392660 

28.23 

3.77 [2.94-6.14] 

0 

0.427 

7 

0 

0 

0 


40.2 J 

3.588210 

-30.392550 

29.01 



0.850 

7 

0 


0 


40.3 J 

3.577540 

-30.399370 

28.53 



0.014 

5 

0 


0 


41.1 J 

3.599170 

-30.399580 

28.45 

4.48 [2.24-6.23] 

1 

1.133 

4 

0 

1 

1 

3.36 [3.24-3.79] 

41.2 J 

3.593790 

-30.407820 

29.78 



1.485 

4 

1 


0 

1.48 [1.42-1.60] 

41.3 J 

3.583460 

-30.408500 

28.40 



0.437 

4 

0 


1 

5.73 [4.78-7.22] 

41.4 J 

3.590420 

-30.404040 

28.90 



0.270 

4 

1 


0 

2.07 [1.75-2.07] 

42.1 J 

3.597290 

-30.400610 

28.56 

3.58 [1.31-6.51] 

1 

1.162 

5 

0 

1 

1 

4.64 [4.35-5.92] 

42.2 J 

3.590960 

-30.403260 

28.69 



0.796 

5 

1 


0 

1.64 [1.55-1.66] 

42.3 J 

3.581580 

-30.408630 

28.54 



0.359 

5 

0 


1 

3.51 [3.59-5.39] 

42.4 J 

3.594250 

-30.406390 

28.58 



3.629 

5 

1 


0 

1.72 [1.67-2.03] 

43.1 J 

3.597830 

-30.402500 

27.35 

3.72 [1.68-6.98] 

0 

0.156 

5 

0 

0 

0 


43.2 J 

3.583960 

-30.409820 

27.22 



0.147 

5 

0 


0 


44.1 J 

3.583460 

-30.406970 

27.75 

0.24 

0 

0.069 

4 

0 

0 

0 


44.2 J 

3.596710 

-30.399760 

27.96 



0.090 

4 

0 


0 


45.1 J 

3.584830 

-30.398470 

28.12 

2.99 [1.77-9.07] 

0 

0.311 

7 

1 

0 

0 


45.2 J 

3.581420 

-30.403950 

27.09 



2.834 

6 

0 


0 


45.3 J 

3.586880 

-30.401290 

26.34 



3.243 

5 

1 


0 


45.4 J 

3.597420 

-30.396150 





4 

1 


0 


46.1 Z 

3.595040 

-30.400750 

27.53 

10.00 6 8.76 [3.40-9.30] 


0.780 

7 

0 

1 

1 

7.58 [6.30-8.67] 

46.2 Z 

3.592500 

-30.401490 

27.89 



0.059 

7 

0 


1 

15.03 [9.21-15.25] 

46.3 Z 

3.577500 

-30.408700 

29.68 



1.565 

6 

0 


1 

2.68 [2.67-2.71] 

47.1 J 

3.586460 

-30.392130 

29.38 

3.58 [1.96-8.93] 

0 

1.843 

2 

0 

0 

0 


47.2 J 

3.585830 

-30.392240 

27.39 



1.952 

2 

0 


0 


47.3 J 

3.589670 

-30.392140 

28.90 



0.423 

1 

0 


0 


47.4 J 

3.578330 

-30.398130 

29.45 



3.755 

1 

0 


0 
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TABLE l:Multiply lensed 

arc 

systems identified 

in 

the Abell 2744 field 




48.1 J 

3.594250 

-30.402840 

27.06 


1.50 [1.45-8.37] 

0 

0.320 

4 

0 

0 

0 


48.2 J 

3.592790 

-30.403130 

28.14 




0.969 

4 

1 


0 


48.3 J 

3.582080 

-30.408580 






4 

0 


0 


49.1 J 

3.592670 

-30.408270 

29.25 


0.04 

0 

0.637 

4 

1 

0 

0 


49.2 J 

3.590210 

-30.408810 






5 

0 


0 


49.3 J 

3.597500 

-30.403140 

28.24 




0.101 

5 

0 


0 


50.1 J 

3.577960 

-30.401620 

28.23 


4.59 [1.12-5.91] 

1 

0.422 

6 

0 

1 

1 

4.01 [3.61-4.09] 

50.2 J 

3.593960 

-30.394300 

28.14 




0.440 

6 

0 


1 

4.07 [4.18-5.19] 

50.3 J 

3.585330 

-30.393630 

26.11 




0.674 

4 

1 


0 

1.93 [1.83-2.05] 

51.1 J 

3.586830 

-30.405550 

28.70 


4.88 [1.91-6.61] 

1 

0.364 

7 

0 

1 

1 

50.72 [26.92-369.66] 

51.2 J 

3.586460 

-30.405670 

29.04 




0.462 

6 

0 


1 

170.12 [35.98-623.76] 

51.3 J 

3.599000 

-30.398300 






4 

1 


0 

2.96 [2.73-3.51] 

52.1 J 

3.586580 

-30.397010 

25.30 


0.01 


0.002 

6 

1 

0 

0 


52.2 J 

3.586130 

-30.397130 






6 

1 


0 


52.3 J 

3.588460 

-30.396860 

29.05 




0.689 

3 

1 


0 


53.1 J 

3.579830 

-30.401590 

28.35 


1.60 [1.64-8.46] 

0 

0.000 

7 

0 

0 

0 


53.2 J 

3.583540 

-30.396700 






7 

1 


0 


53.3 J 

3.597040 

-30.394550 






1 

0 


0 


54.1 J 

3.592370 

-30.409890 

27.59 


6.69 [1.59-8.56] 

0 

0.294 

7 

0 

0 

0 


54.2 J 

3.588250 

-30.410340 






7 

1 


0 


54.3 J 

3.588420 

-30.410320 

27.25 




0.289 

7 

1 


0 


54.4 J 

3.590080 

-30.410270 






7 

1 


0 


55.1 L 

3.597042 

-30.404746 

25.57 

1.50 

0.03 


1.777 

8 

0 

1 

1 

12.63 [10.24-24.00] 

55.2 L 

3.596371 

-30.406161 

25.84 




1.463 

8 

1 


1 

10.28 [4.69-9.87] 

55.3 L 

3.585726 

-30.410085 

27.14 




1.417 

8 

0 


1 

3.50 [3.07-3.45] 

56.1 G 

3.582526 

-30.402290 

25.52 

1.20 

2.61 [1.08-6.09] 


1.934 

8 

0 

1 

1 

5.54 [4.09-7.62] 

56.2 L 

3.596728 

-30.396298 

26.88 




2.895 

8 

0 


1 

2.52 [2.53-2.61] 

56.3 L 

3.584467 

-30.399292 

26.52 




4.941 

8 

1 


1 

10.39 [9.14-11.08] 

56.4 L 

3.586219 

-30.400850 

24.28 




4.457 

8 

1 


1 

4.04 [3.57-4.56] 

57.1 I 

3.598676 

-30.40491 

24.59 


1.05 [0.71-2.70] 

0 

0.490 

0 

0 

0 

0 


57.2 I 

3.587053 

-30.41126 

24.98 




0.733 

0 

0 


0 


57.3 I 

3.596818 

-30.40783 

27.40 




0.494 

0 

1 


0 


58.1 I 

3.578090 

-30.39964 

25.04 


0.01 

0 

7.221 

4 

0 

0 

0 


58.2 I 

3.589237 

-30.39444 

24.18 




3.590 

4 

0 


0 


59.1 I 

3.584284 

-30.40893 

26.84 


4.27 [1.45-6.56] 

1 

0.037 

1 

0 

1 

1 

8.12 [4.82-9.31] 

59.2 I 

3.598125 

-30.40098 

27.61 




0.108 

1 

0 


1 

5.15 [3.47-6.64] 

60.1 G 

3.5980768 

-30.403991 

27.39 


0.07 

0 

2.501 

8 

0 

0 

0 


60.2 G 

3.5957195 

-30.40755 

27.27 




0.763 

8 

1 


0 


60.3 G 

3.587391 

-30.410149 

27.37 




1.588 

8 

0 


0 



Note. — Objects for which F140W magnitudes are not reported are not detected by SExtractor. Contaminated objects are only required to fulfill 
the morphology criterion, M > 5, but they are not used to compute ^Bayes- Arc images with pass flag =§ 1 fulfill Color+Morphology+Contamination 
criteria & ^Bayes > ^cluster & ^Bayes single-valued. Note that arc systems 1 and 55 have the same physical origin and therefore have the same 
2 sp ec through our identification of arc 1.3 (object ID #336). So do systems 2 and 56, for which z spe c is measured through arc 56.1 (object ID 
#888). Systems 15, 16, and 17 are not included in the table or the lens model because they belong to northern subclumps with > 1 arcminute 
separation from the cluster center shown in Fig. 3. The coordinates of these systems can be found in Richard et al. (2014). 

1 G = this work, J = (Jauzac et al. 2014), T = (Johnson et al. 2014), R = (Richard et al. 2014), C = (Clement et al. in preparation), Z = (Zitrin 
et al. 2014), L = (Lam et al. 2014), I = (Ishigaki et al. 2015). References correspond to most recent quote in the literature. System 60 is identified 
in this work for the first time. 

2 

Redshift obtained from hierarchical Bayesian modeling. The values outside of the brackets are the mode of the combined posterior probability 


distribution. For systems that fulfill the selection criteria and do not have a spectroscopic redshift, this is the redshift assigned to the system in 
the lens model. Uncertainties represent the 68% credible region. Note that ZBayes = 0-01 is assigned if the posterior distribution of photometric 

redshift declines monotonically from 0, and is thus considered highly uncertain. 

3 

Quality 0 indicates that ZBayes is unreliable due to ^Bayes < ^cluster and/or there exists strong multi-modality in the posterior probability 
distribution of photometric redshift and/or only one image was used to compute ZBayes- Quality 1 indicates that ZBayes is secure. 

^ Best-fit magnification and 68% confidence limits derived from resampling the multiple image systems themselves and their photometric redshifts 


from the combined posterior distributions. 

5 

Due to the use of a fixed SExtractor detection image at F160W, 2.2 was not detected with even the most aggressive SExtractor detection 


settings, i.e. the “hot” mode settings. Upon visual inspection in other HST bands, however, the object is clearly separated and unmistakably 


belongs to system 2. 

^ The redshift of this system comes from the geometric constraint by Zitrin et al. (2014). 












TABLE 2:Emission line detection results on multiply and singly lensed sources 
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Note. — The first part of this table consists of emission line identifications for the arcs of quality levels 4, 3, and 2, whereas the second part is comprised of only high-confidence (quality 2 
4) emission line objects newly discovered during the blind search procedure, as described in Section 4.3. The uncertainty ranges and error bars quoted here all represent 1-cr confidence level. 

1 Magnifications of multiply-lensed objects are calculated assuming redshift z spec , which was only used in the lens model for quality> 2 objects. 
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Fig. 16.— : Left: System 60, an example of a multiple image system whose images fulfill the morphology criteria, all 
three of which receive maximum grades (M) of 8. The images all possess the same point-like morphology and have 
similar surface brightnesses. This multiple image system is discovered for the first time in this work. Unfortunately, 
the photometric redshifts for this system were not reliable, so we could not use the image system in the lens model. The 
RGB image for system 60 is made from a combination of frames: F105W, F125W and F160W, with a pixel scale of 60 
milli-arcseconds. Right: System 27, a multiple image system whose images fail to fulfill the color/morphology criteria 
due to poor morphology grades. All three images receive a morphology grade of 5. Both inspectors independently gave 
the images a low grade because they were too faint to discern any shareable features. The color scale is the same for 
all three images, showing the severe foreground contamination due to the ICL in image 27.3. Without spectroscopic 
redshifts and morphological details, there is no assurance that these images are of the same source other than their 
colors, which are largely uncertain due to the extreme faintness of the images. The RGB image for system 27 is 
made from a combination of frames: F435W, F814W and F125W, with a pixel scale of 60 milli-arcseconds. Out of all 
possible RGB combinations of HST images, the images are most discernible in this one. In both system 60 and system 
27, the radius of the green circles is 0.72 arcseconds. 



Fig. 17.— : Left: System 7, an example of a multiple image system whose images receive high morphology grades (M) 
of 8, 7 and 8, respectively. The RGB image for system 7 is made from a combination of frames: F105W, F125W and 
F160W. Right: System 57, a multiple image system that we can rule out due to an impossible lensing configuration, 
for which M = 0 is assigned for all images. The proposed images are in a naked-cusp configuration and should 
all therefore have similar magnification. 57.3 is barely detected, yet 57.1 and 57.2 are detected at high S/N, thus 
invalidating the system. If 57.3 were simply misidentified, there should still be a counter-image with the same color, 
similar magnification and elongation as 57.1 and 57.2 nearby 57.3. We find no such counter-image. There may be 
another lensing configuration in which 57.1 and 57.2 are indeed multiple images of the same source, in which case 
there may be additional counter-images elsewhere in the cluster field that have yet to be found. The RGB image for 
system 57 is made from a combination of frames: F435W, F606W and F105W. In both systems 7 and 57, the radius 
of the green circles is 0.72 arcseconds. 
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Fig. 18.— : Two examples of arc systems that we exclude from the lens model due to poorly-constrained redshifts. 
Shown are the photometric redshift probability density functions, P(z), for the individual images (blue, green, red) 
along with a combined P(z) obtained via hierarchical Bayesian modeling (purple). Left System 38 is excluded because 
it is bimodal with comparable probabilities under each peak. Right System 40 is excluded because the individual 
photometric redshifts of images 40.1, 40.2 and 40.3 are in significant disagreement. 
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Fig. 19.— : Two examples of arc systems that we exclude from the lens model due to a peak in the combined Bayesian 
redshift below the cluster redshift of Abell 2744, z = 0.308. Left System 9. The P(z) for image 9.1 is not shown 
because the object is flagged as contaminated. Right System 20. Likewise, P(z) for image 20.2 is not shown because 
the object is flagged as contaminated. 



































